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We consider wires near a zero temperature transition between superconducting and metallic states. 
The critical theory obeys hyperscaling, which leads to a universal frequency, temperature, and length 
dependence of the conductance; quantum and thermal phase slips are contained within this critical 
theory. Normal (NN), superconducting (SS) and mixed (SN) leads on the wire determine distinct 
universality classes. For the SN case, wires near the critical point have a universal d.c. conductance 
which is independent of the length of the wire at low temperatures. 



The fluctuations of superconducting order in wires 
have long been the focus of experimental interest, and 
recent measurements Q, |^ have extended such obser- 
vations to the nanoscale. Such wires have a diameter 
which is of the order of the BCS coherence length or 
larger, so that there are a large number of transverse 
channels for electronic conduction and the single elec- 
tron levels are effectively three-dimensional. However, at 
low energies the collective fluctuations of the supercon- 
ducting order are one-dimensional because the diameter 
of the wire is much smaller than its length, L. Above, 
and not too far below the bulk superconducting critical 
temperature Tc, these measurements have been success- 
fully interpreted using a theory ^ |^ based upon the 
time-dependent Ginzburg-Landau (TDGL) equation. At 
very low temperatures (T), there is a crossover from the 
purely classical thermal fluctuations of TDGL theory, to 
effects associated with quantum fluctuations of the su- 
perconducting order. In particular, as the normal state 
resistance of the wires is increased, they apparently un- 
dergo a superconductor to metal quantum transition. In 
superconducting wires, there is superflow with infinite 
conductance as T 0. In contrast, metallic wires have a 
finite conductance, g, as T —f which decreases inversely 
with L, g = a / L, where a is the conductivity. 

This paper will present new results on the transport 
properties of wires in the vicinity of the T = quantum 
superconductor-metal transition (SMT). We will argue 
that the conductance has a singular contribution which 
is a universal function of L, T, and the measurement fre- 
quency, iu, as specified in Eq. ^ below. At a formal level, 
this scaling form parallels that proposed earlier for the 
superconductor-insulator transition (SIT) in d = 2 spa- 
tial dimensions thermodynamic limit; we con- 
tend here that such scaling arguments can be extended 
to the SMT in d = 1, and to L finite (even though they 
cannot be extended to the SIT in d = 1). Furthermore, 
for L finite, we use the theory of surface critical behavior 
to demonstrate that the leads connected to the sample 
determine distinct universality classes of the conductance 
scaling function: wires with superconducting (SS), nor- 
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FIG. 1: A wire which is tuned from a superconductor to a 
metal by (say) reducing its diameter. The leads on the wire 
are either normal (N) or superconducting (S) and the NN, 
SN, and SS cases belong to distinct universality classes. 

mal (NN), and mixed (SN) leads have distinct scaling 
functions, but other details of the leads arc unimportant 
(see FigCJ. For wires close to the SMT with L < Lt, 
where Lt ~ T~^/^ is a thermal cutoff length specified 
later (z is the dynamic critical exponent), our arguments 
imply that as w ^ 0, .g = 'ie'^Cx /h (with X = NN, SN, 
or SS), where the Cx are universal numbers. Thus the 
d.c. conductance of the wire is independent of L — this 
happens because the physics is dominated by low energy 
superconducting fluctuations whose characteristic size is 
L itself. We will determine the values of Cx in a large 
n limit (n is the number of real components of the order 
parameter, and the case of interest here is n = 2), and 
by quantum Monte Carlo simulations for a non-random 
universality class. 

The framework of our analysis is provided by a theory 
for the T = SMT proposed by Feigel'man and Larkin 
[i(tI | , and examined in a number of studies [HI in 
d = 2. This theory may be viewed as a natural quantum 
extension of the TDGL — it reduces to the TDGL at high 
T. We suspect it is also the theory of critical fluctuations 
in the analysis of suppression of the critical temperature 
m RefJlJ The key to our conclusions is the demonstra- 
tion |lfil | that this theory obeys conventional hyperscaling 
properties at the T = SMT in d < 2 in the non-random 
class. Transport properties in the vicinity of the d = 1 
SMT are controlled entirely by the critical theory, and 
perturbations from irrelevant operators are not needed. 
In this respect, the situation is similar to the SIT in d = 2 
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[g| . It is important to note that scaling forms like Eq. 
do not apply to the non-random SIT in d = 1 because 
this SIT is in the Kosterlitz-Thouless universality class, 
and the conductivity of the superfluid phase near the SIT 
is determined by 'irrelevant' phase-slvp operators |l7lll8| . 

We now state our central scaling hypothesis for the 
d = 1 SMT. The conductance always has a 'background' 
contribution from a parallel metallic channel, and so we 
write g = (Ji,/ L + 'g, where ai, is the background conduc- 
tivity which is not expected to have any important L, 
T, or Lo dependence; note that the background contribu- 
tion to the conductance always falls as 1/L, and so can 
become negligible compared to the singular contribution 
g. The latter contribution is universal, and obeys the 
scaling form 



described by the imaginary time (r) partition function 
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Here $x is a universal scaling function (note that the 
overall scale of is universal and there is no non- 
universal prefactor), v is the correlation length expo- 
nent, w is the parameter which tunes the wire (say, its 
diameter) across the superconductor-metal transition at 
w — Wc, and ci,2 are (the only) non-universal constants. 
For w > Wc Eq. describes metallic conduction, and 
for w < Wc quantum and thermal phase slips disrupt 
the superflow; we emphasize that, unlike the Kosterlitz- 
Thouless SIT theory [l7|, such phase slips are contained 
within the critical theory of the SMT. In this first dis- 
cussion here, we focus on the quantum critical behav- 
ior of short wires with L < Lt (cifcsT)^^/^ and 
L < {c2\w — Wc\)^'^ ■ In this case we can write Eq. ^ as 



giiu) = iAeyh)Fxiy) 



(2) 



where Fx{y) = ^x{y, 0, 0). Computations of the univer- 
sal function Fx {y) are provided below. 

We orient ourselves, and estimate various parameters, 
by recalling the TDGL approach, following the notation 
of Ref. 0. The spatial {x) and temporal (t) evolution of 
the complex superconducting order parameter 'i'{x,t) is 
determined by the classical equation of motion 



hjdt-^ = -[a + - Sdl] * + / 



(3) 



where / is a Langevin noise term, a = ao{T — Tc)/Tc, and 
ao, 6, (5, 7 are T-independent parameters whose values are 
specified in Ref. 0. The dissipative co-efficient 7 arises 
from the decay of Cooper pairs into normal electrons. 
The considerations of Refs. llOllllI show how a quantized 
version of Eq. (0) can also apply near a T = SMT in 
systems with an inhomogeneous BCS coupling between 
the electrons, with a reservoir of normal electrons being 
provided by regions of the sample with a weaker bare 
coupling (here, these could be near the wire surface). 
For a wire in the region < x < L, such a theory is 
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Q|*(0,r)p + a|*(i,T)|2 



(4) 



-Rc[i7£«'(0,T)] -Re[i/^^'(L,r)] . (5) 



Here A is the cross-sectional area of the wire, (3 = 
h/{kBT), ujn is a bosonic Matsubara frequency, and a 
tunes the system across the SMT, and so a = ao{wc — 
w)/wc- As a first guess, we may estimate that the param- 
eters ao, b, S, and 7 have the same values as those esti- 
mated in Ref. in the dirty limit, although the value of 7 
will decrease at low T due to reduced damping. The pres- 
ence of disorder in the wire also implies a quenched ran- 
dom X dependence of all the couplings in 5buik; our quan- 
titative results below are limited to the non-random uni- 
versality class where such x dependence is neglected. The 
term is permitted in .Sbuik but its co-efficient is 

proportional to the degree of particle-hole asymmetry on 
the scaling of pairing energy, and should be quite small: 
we defer analysis of its consequences to later work. The 
terms in iSboundary reflect the presence of the left /right 
(l/r) leads, with C encodin g th e boundary conditions on 
the superconducting order [igj. For a N lead we have 
H = and C > 0, while a S lead has H Q because the 
bulk superconductivity of the lead acts like an boundary 
ordering field on ^E* in the wire. 

We now discuss the properties of the QMT of iSbuik 
in the thermodynamic limit. These were described re- 
cently by Pankov et al. The QMT has an up- 
per critical dimension of d = 2, and universal critical 
properties were computed in an expansion in e = 2 — d. 
This expansion obeys hyperscaling properties to all or- 
ders in (2 — d), and justifies the scaling assumptions be- 
hind Eq. The long-range interaction between 
^I* fluctuations generated by the |a;| dissipative term is 
preserved under renormalization, and this leads [isl 
to the exponent identity z — 2 — rj, where ry is the 
anomalous dimension of ^> which was computed to be 
ri = (n-l- 2)(12 - n'^)e^/{A{n + 8)'^) to order e^. We com- 
puted v by similar methods to the same order and ob- 
tained 
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This computation also provides the mean field estimate 
Lt ~ a/ 5/ {'-/kBT). We have also carried out quantum 
Monte Carlo simulations on a lattice realization of iSbuik 
(described below), following those of the n = 1 case in 
Ref. (23. The results for the exponents are similar, with 
z = 1.97(3), z + 77 = 1.985(20) and v = 0.689(6). These 
are in good agreement with the predictions of the e ex- 
pansion which upon evaluation for n = 2, e = 1 yield 
rj = 0.021 and v = 0.663 (from the scries for or 

= 0.701 (from the series for l/i^). 

The influence of 5boundary can be understood using the 
theory of surface critical behavior 9]. A N lead cor- 
responds to the ^ordinary transition': C is a relevant 
perturbation which flows to C = cx) under renormaliza- 
tion, and so we have Dirichlet boundary conditions with 
^'(0; L, r) = in the scaling limit. Similarly a S lead cor- 
responds to the ''extraordinary transition' in which the 
magnitude of the ordering field \H\ scales to 00. In both 
cases, the structure of 4* correlations near the edge is uni- 
versal, and independent of the specific values of H and 
C. For the SS case there will be a residual universal de- 
pendence on the phase difference A<i> = a.rg{H^ Hr); we 
take A<i> = in our computations below. 

We have now assembled the tools needed to compute 
the conductance using the Kubo formula '2l!| . We define 



K{ujn) = / dT (J(t) J(0)) e'' 
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where J(t) = (SA/itiLi)) dx{'<i!*d^'i' - 9^***) 
and D ^ {2dA/{hL^))J^dx{\'f\^). Then g{uj) = 

A first computation of the conductance was obtained 
in the n = 00 limit. The bulk theory has the expo- 
nents z = 2, rj ~ V = 1. The saddle point equations 
I22II23I I for finite L involve determination of the optimum 
a;-dependent values of the decoupling field for the quar- 
tic term at the bulk quantum critical point, and this was 
done numerically after discretizing x to a lattice. The 
saddle point solution was inserted into Eq. ITJ, and leads 
to the results shown in Fig|2| Note the remarkable insen- 
sitivity of the conductance to large changes in the values 
of the boundary couplings, which confirms our analysis 
of the surface critical behavior for large n. Universality 
is also demonstrated in Fig O where distinct couplings 
scale to the same conductance in the large L limit. 

These results show interesting differences in the small 
Lo behavior of the conductivity for the NN, SN and SS 
cases. In the NN case we find Re[JNN(y)] = i^iv"^ + ■ ■ ■ 
for small y (the k co-efficients here and below are univer- 
sal numbers), and hence the universal d.c. conductance 
vanishes with Cnn = 0. This is a consequence of the 
limited density of states for decay of Cooper pairs at low 
energies, and we expect this feature to hold beyond the 
large n expansion (it is also consistent with our Monte 



0.20 

0.15 
0.10 



X=SN 



//■ 

;/■/ 
/V 

V 



(a) 




FIG. 2: The universal function Re [Fx] in the large n limit 
for real frequencies, with y = ^hwh"^ /5. We discretized x in 
Eq. ||7|l on a lattice of spacing unity, and rescaled r and ^ 
to obtain 5A/Ti = 1 and A7 = 1, and set 6 = 1 and used 
a ultraviolet frequency cutoff tt. (a) Test of universality for 
the SN case with L = 400. The parameters {Ci,Cr,He,Hr) 
are (i) (1,1,1,0), [it) (1,1,10,0), (m) (1,10,1,0), (iv) (5,1,5,0), 
(v) (1,1,0.5,0). A similar insensitivity to boundary param- 
eters was found for the NN and SN cases, (b) Results for 
all 3 universality classes are obtained for L — 800, other 
parameters as in Figl^Ji, and {Ce,Cr, Hi, Hr) taking values 
(1,1,0,0) for NN, (1,1,1,0) for SN, and (1,1,1,1) for SS. For 
the SS case, there is an additional Josephson current con- 
tribution Fss(y) = TrQS{y), with g — 2.98, which is not 
shown above. All three classes have the common behavior 
Fx{y ^ CX)) = 0.81994(1 + i)/^ for large n. 



Carlo results below). In contrast, for the SN case we find 
Re[i^SN(2/)] = CsN + K2?/^ -I- . . ., with Csn non-zero and 
universal. Here the proximity effect of the S lead provides 
a new channel for metallic conduction. Finally, for the 
SS case we obtain Re[Fss(2/)] = '^Q^{y)+Css + ^^sV^ + . . .. 
Here there is a residual Josephson coupling, proportional 
to the universal number g, between the S leads, induced 
by the proximity effect on both leads. 

Finally, we describe our quantum Monte Carlo sim- 
ulations on a "hard spin" lattice realization of Z. We 
discretize space into points j = 1...L (integer), imag- 
inary time into points k = \ . . . [3 (integer measuring 
fi/{kBT)), and set ^'(xj, t^) = e"'^^-'' to a unit modulus 
complex number. The continuum theory Z is realized by 
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FIG. 3: Extrapolation of the d.c. conductance Csn to the uni- 
versal scaling limit (L oo). The large n parameters are as in 
FigUl but with {b,Ct,Cr,Ht,Hr) = (i) (46.25,1,1,0.745,0), 
(m) (0.925, 1, 1, 7.45, 0). The quantum Monte Carlo parame- 
ters are as in the text, with {He,Hr) = (1) (10,0), (2) (1,0). 

^0 = Ilj=i nLi lo'' '^^3,k exp(-50) with 

L-l [3 

S4, = -~Kx ^ ^ COS {4>j^k - <t>] + l,k) 
J = l fc=l 
L P 

- X! X! ^^i^^ ~ '^i.fc+i) 

j=l k=l 

P 

- ^ [Hi COs{(f>i^k) + Hr COs{(f>L,k)] 
fc=l 

I " /^^ V [1 - COS (03,fc - '/>j,fc')] 

Bin^Wn.-r.0//3)] ' ^ ^ 

where temporal periodic boundary conditions are im- 
posed by identifying (/'j,/3+i = 4>i,i- The Kr term be- 
comes l^r^'P in the continuum limit: such a term is 
formally irrelevant and so was not explicitly noted in 
Z. Note that the a damping term in derives from 
the 7 term in 5buik, and it differs from that usually 
assumed in the phenomenological resistively-shunted- 
Josephson (RSJ) junction models 0) Ell — it depends 
only on cosines of phase differences, while that in the 
RSJ models depends upon squares of phase differences; 
this feature is crucial to our results. We chose a — 0.3, 
Kr = 0.1, and determined the bulk critical point to 
be at Kx = 0.92132(2). These values were used in 
the subsequent computation of the conductance from 
Eq. O with J(rfc) = {Kx/L) Yl!jZl iiva(4>j.k-4>3+i,k) and 
D = (Kx/L^) X;t=i^ (cos(0j, fc - (f>j+i,k))- 

We obtained Monte Carlo results for Fx{y) along the 
imaginary frequency axis, and the results had a structure 
similar to those of the large n theory. We show in Fig El 
the values of the universal d.c. conductance Csn as a 
function of 1/i: the large n theory is seen to significantly 



underestimate its value, but has a similar sensitivity to 
finite sizes. 

This paper has described the consequences of a the- 
ory 0, O, E, for a quantum transition between 
a superconductor and a metal in one spatial dimension. 
Our results apply to wires in which the superconducting 
'Cooperon' fiuctuations are effectively one-dimensional, 
but there are a very large number of transverse single- 
electron channels so that the strictly one-dimensional 
Luttinger liquid physics does not apply. We have argued 
that the proliferation of thermal and quantum phase slips 
near such a transition is conveniently described by a 'soft- 
spin' continuum theory in Eqs. (|4I5|I . whose critical point 
obeys conventional hyperscaling properties. We used an- 
alytic and Monte Carlo computations to make predic- 
tions on a universal conductance. We hope that future 
experiments on short wires can test our predictions, and 
partic ularly _ their sensitivity to the leads; other recent 
works 25, 26] have also addressed lead-sensitivity, but in 
very different frameworks. 
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